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Abstract 

The logarithmic and constant contributions to the Wilson coefficient of the longitudinal 
heavy quark structure function to O(a^) are calculated using mass factorization techniques 
in Mellin space. The small x behaviour of the Wilson coefficient is determined. Numerical 
illustrations are presented. 
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1 Introduction 



Deeply inelastic electron-nucleon scattering at large momentum transfer provides one of the 
cleanest possibilities to test the predictions of Quantum Chromodynamics (QCD). In the case 
of pure photon exchange the structure functions F 2 (x, Q 2 ) and F L (x, Q 2 ) describe the scattering 
cross section. While the former structure function is well measured in a wide kinematic region [1] , 
F L (x,Q 2 ) was mainly measured in fixed target experiments [2] and determined in the high 
y region at HERA [3] using an extrapolation method. Future detailed measurements of the 
longitudinal structure function Fl(x,Q 2 ) at HERA are still to be performed [4]. At leading 
order in the coupling constant the gluon distribution g(x, Q 2 ) does not contribute to the structure 
function F 2 (x, Q 2 ) directly, but only to its derivative, which weakens the sensitivity. In the region 
of smaller values of x the structure function F L (x,Q 2 ) = F 2 (x,Q 2 ) — 2xFi(x,Q 2 ), however, is 
dominated by the gluon contribution. Therefore, this structure function may yield essential 
constraints on g(x, Q 2 ). In lowest order in the coupling constant (a° s ) and vanishing target-mass 
effects, the twist-2 contributions to the structure functions F 2 and F L obey the Callan-Gross [5] 
relation 

F 2 (x,Q 2 ) = 2xF 1 (x,Q 2 ), F L (x,Q 2 ) = . (1) 

F L (x,Q 2 ) receives leading order contributions due to target mass effects [6]. The Callan-Gross 
relation is further broken by QCD corrections. The corresponding Wilson coefficients for massless 
quarks were calculated in leading (LO) [7], next-to-leading (NLO) [8-10], and next-to-next-to- 
leading order (NNLO) [11-13]. Since the leading order coefficient functions are polynomial, 
scheme-invariant quantities one may construct a simple mapping of F^°(x,Q 2 ) to g LO (x, Q 2 ) 
taking the quark distributions from the F 2 measurement [14]. The leading small x terms for 
the coefficient functions of Fl have been derived in [15] and agree with the known fixed order 
results (NLO, NNLO) [8,9,12,13]. The gluonic contribution to F L (x : Q 2 ) was calculated using 
the k± representation in leading order [16], which turns out to be numerical very close to the 
NLO result [8]. The numerical impact of the small x resummation [15] on Fl was studied in [17]. 
Similar to the small x resummation for the splitting functions, formally sub-leading terms lead 
to comparable but widely compensating effects, as seen comparing the magnitude of these terms 
for fixed orders in the coupling constant. This behaviour was later observed also in [13]. To 
draw firm conclusions on the effect of these resummations, several sub-leading series of terms 
have to be known. Higher twist contributions to F L (x,Q 2 ), partly under model assumptions, 
were considered in [18]. 

Since the longitudinal structure functions Fl(x,Q 2 ) contains rather large heavy flavor con- 
tributions in the small x region [19], a consistent analysis has to account for these effects, which 
were calculated in leading [20] and next-to-leading order [21, 22]. 2 The NLO corrections [21] 
could not be performed in analytic form completely. This is also expected for even higher or- 
ders, due to the complexity of the phase space integrals. However, complete analytic results 
may be derived in the asymptotic region Q 2 ^> m 2 calculating all contributions but the power 
suppressed terms (m 2 /Q 2 ) k , [24, 25]. 3 

In the present paper we use the method of Ref. [24] to derive the heavy quark Wilson 
coefficients for F® Q (x,Q 2 ) to 0(a^) in the region Q 2 ^> m 2 . In Section 2 we give a brief 
outline of the method. The Wilson coefficients are derived in Section 3. Their small x behaviour 

2 Fast Mellin-space expressions for these Wilson coefficients were given in [23]. 
3 For related work for other processes, see [26]. 
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is discussed in Section 4. In Section 5 numerical are presented and Section 6 contains the 
conclusions. Some useful relations are summarized in the Appendix. 



2 The Method 

In the twist-2 approximation the nucleon structure functions Fi(x,Q 2 ) are described as Mellin 
convolutions between the parton densities fj(x,fj, 2 ) and the Wilson coefficients Cj(x, Q 2 / n 2 ) 

F t (x, Q 2 ) = £ C? (x, ® fox, /i 2 ) (2) 



to all orders in perturbation theory due to the factorization theorem. Here /i 2 denotes the 
factorization scale and the Mellin convolution is given by the integral 

[A <g> B] (x) = [ dx ± [ dx 2 S(x -x ± x 2 ) A(x 1 )B(x 2 ) ■ (3) 

JO JO 

Since the distributions fj refer to massless partons, the heavy flavor effects are contained in the 
Wilson coefficients only. We are interested in the massive contributions in the region Q 2 ^> m 2 . 
These are the non-power corrections in m 2 / Q 2 , i.e. all logarithmic contributions and the constant 
term. We apply the collinear parton model, i.e. the parton 4-momentum is p = zP, with P 
the nucleon momentum. The massive Wilson coefficients itself can be viewed as a quasi cross 
section in pl / *-scattering, where V* denotes the exchanged virtual vector boson. In the limit 
Q 2 ^> m 2 the massive Wilson coefficients H 2 j* ti (Q 2 /m 2 ,m 2 j /i 2 ,x), likewise the case for the 
structure functions (2), factorize 

"is(t^A=cs(t*)®<r(^A (4) 



into Wilson coefficients C (Q 2 / /u, 2 , x) accounting for light flavors only and massive opera- 
tor matrix elements A^f s (m 2 /fi 2 , x). The latter take a similar role as the parton densities 
in (2), but are perturbatively calculable. The factorization (4) is a consequence of the renor- 
malization group and the fact that we restrict the investigation to non-power corrections to 
H L ^ (Q 2 /m 2 , m 2 /yU 2 , x). The operator matrix elements A^ s obey the expansion 

A S,NS ( m2 \ /•| /n S,NS|-\ r . / ^S,NS,(0 • /rN 

\ ' / 1=1 

S NS 

of the twist-2 quark singlet and non-singlet operators O k ' between partonic states which 
are related by collinear factorization to the initial-state nucleon states \N) and a s = a s (/z 2 )/(47r) 
denotes the strong coupling constant. The operator matrix elements are process-independent 

S NS 

quantities. The process dependence of H L \ is described by the associated coefficient functions 



^ff)=E^(fV fc = NS,S,s. (6) 



The MS coefficient functions, in the massless limit, corresponding to the heavy quarks only, are 
denoted by 

= C L , k N L + N H ) - C Lik N L ) , (7) 
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where Nh,Nl are the number of heavy and light flavors, respectively. In the following we will 
consider the case of a single heavy quark, i.e. N H = 1. The formalism is easily generalized 
to more than one heavy quark species. The heavy flavor Wilson coefficient is obtained as the 
expansion of the product of (5,6) to the respective order in a s . 



3 The Wilson Coefficients in the Region Q 2 > M 2 

In the limit of vanishing nucleon mass effects, cf. [6], the longitudinal structure function emerges 
only at 0(a s ) due to the Callan-Gross relation [5]. The leading order contribution is purely 
gluonic [71. At 0(a 2 ) F? Q (x,Q 2 ) receives also quarkonic contributions. 

' r S,PS,NS 



To 0(a 3 s ) the heavy quark Wilson coefficients read 
f Q 2 m 2 \ 



rrPS 

lj,q 



ttNS 

tlj n 



where 



9 ' 9 



'Q 2 

m 2 ' fi 2 

'q 2 rrf 

m 2 ' fx 2 



A 



(i) 

Q,9 { m 2 



+ < 



A (2) ifl' 
A Q,9 m 2 



Q,g 



2 ^PS,(2) / Q 2N 

2^NS,(2) / Q 2 ^ 



A PS,(2) / V_ 



,NS,(2) / }J_ 

gg,Q \ m 2 



C 



(2) 
L,q 



PS 
L,q 



m 



/i 



, r (l) (Q 2 \ , ANS,(3) / <5 2 " 



(8) 
(9) 

(10) 
(11) 



and 



lj,q 



rrNS , rrPS 
H L,q + ^L, q ■ 



(12) 



C^J(Q 2 1 /i 2 ) are the scale dependent Wilson coefficients in the MS scheme with C^J(Q 2 /fi 2 ) 

c ( l\ for Q 2 = fi 2 given in [7-13]. 

The operator matrix elements were derived in [24] and read : 

,2' 



A 



(i) 
Qg 



2 qg U 2 , 



l Qg 



41 = i{P™®[P™-P™ + 2 i3 „]}l„ 2 



m 



(13) 



4 PS,(2) 

Qq 



,NS,(2) 



4R 1)+ °SK ) -^° )+2 / 3 o]}ln(^ 



+77 (1) 



^ 0) - ^ + 2ft 



--Pi 0) ®Pl 0) ln 2 



+ a, 



(2) 

Qg 



Qg 



gq 



V A* 

+„ PS '( 2 ) + o« p(o) 
A),o -™ - 9 /m° N 



BPS,(i) _ „(i)p(o) 



In 



p(0) l n 2 f m \ 1 BNS,(1) [ m \ , NS,(2) 1 fl > p0 



// 



(14) 

(15) 
(16) 



with 



/ = f(N F + 1) - f(N F ) . (17) 
For later fast numerical representations we express the above functions fi(x) in Mellin space, 



M[f t (x)](N) 



dx x N ~ l fi(x) 



(18) 



at (even) integers JV and arrange for analytic continuation to complex values of iV starting from 
these values. 

The splitting functions are 



-2S!(N- 1) + 

N 2 + N + 2 
" !l v ' JV(JV + l)(iV + 2) 



(7V-1X37V + 2) 
2JV(JV + 1) 



n ( ?(Ao 



8T R JV F - 



AC F 



-S 1 (N-1)- 
N 2 + N + 2 



JV 3 — 3JV — 4 



(JV-l)JV(JV + l)(JV + 2) 



+ 2A) 



(JV — 1)JV(JV + 1) 

5JV 5 + 32iV 4 + 49 iV 3 + 38JV 2 + 28JV + 8 



ii ^ > UOF R (AT- 1)JV 3 (JV + l) 3 (JV + 2) 2 

(1) W = ^ S ' (1) = C F T R {±f Sl (N-l)- 3 p 2 (N-l) 

4 (JV — 1)(3JV + 2)(JV 2 - 11JV - 6) 
~9" 



8C F T fi 2 



JV 2 (JV+ l) 2 
JV 2 + JV + 2 Fo2 



(19) 
(20) 
(21) 
(22) 
(23) 

(24) 



Si(N)-S 2 (N)\ -— S 1 (N) 



+ 



JV(JV + l)(JV + 2) 
5JV 6 + 15JV 5 + 36JV 4 + 51JV 3 + 25JV 2 + 8JV + 4 



JV 3 (JV + l) 3 (JV + 2) 



+IQC A T R ■ 
+4 



JV 2 + JV + 2 



JV(JV + l)(JV + 2) 
2JV + 3 



S 2 (JV) + S 2 (N) - C 2 - 2/3'(JV + 1) 



(JV + l) 2 (JV + 2) 2 



^i(JV) + 



A(JV) 



(jv — i)jv 3 (jv + i) 3 (jv + 2y 



(25) 



where 



JV 9 + 6JV 8 + 15JV 7 + 25JV 6 + 36JV 5 + 85JV 4 + 128JV 3 + 104JV 2 
+64JV + 16 . 



The expansion coefficient of the /3-function for the case of light and heavy (Q) flavors read 



A) = ^C A -^T R N f , 



Pi 



0,Q 



-Tr- 



(26) 

(27) 
(28) 



The Mellin transforms lead to harmonic sums. Their analytic continuation for single harmonic 
sums is given by 



S 1 (N-1) = ip(N) + 7e 



S k (N-l) = 
S^(N-l) = 
S. k (N-l) = 



(k-iy. 

I 1\N-1 



^k+N 



^ k - l) (N) + ( k , k>2 
P(N) - ln(2) 



(3 {k \N) 



2k- 



Ck, k>2. 



(29) 
(30) 

(31) 
(32) 



Here ip(z) = d\n[T(z)]/dz, ^ E denotes the Mascheroni-Euler number, ( k the values of Riemann's 
^-function for integer k > 2 and 



m = 5 



(33) 



Multiply nested harmonic sums are reduced to Mellin transforms of basic functions [27-29] for 
which the analytic continuation to complex values of N is performed [30,31]. 

The functions emerging in the scale independent contributions of the operator matrix ele- 
ments are 



1 qV 







a { g{N) = AC F T R < 



N 2 + N + 2 



(34) 
(35) 



~Sl(N-l) + p 3 (N-l) 



N(N + 1) (N + 2) 

—Si(N - 1)S 2 (N - 1) - 2C 2 5 1 (iV - 1 
TV 4 + 16 TV 3 + 15 N 2 - 8 N - 4 



+ 



N(N + 1) 



S 2 (N - 1) 



S 2 (N-l) 



+ 



N 2 (TV + l) 2 (TV + 2) 
3 iV 4 + 2/V 3 + 3iV 2 -4iV-4 

2N 2 (TV + l) 2 (TV + 2) 
N 4 — iV 3 — 1QN 2 + 27V + 4 

+ TTT-T , 2/Ar , ^ S 1 {N - 1) + 



<2 



+4C A T R 



TV 2 (iV + 1) (AT + 2) 

( N 2 + N + 2 
\N(N + l)(N + 2) 



2/V 4 (7V+ l) 4 (/V + 2) 





Li 2 (x) 


AM 






1 +x 



(N) + -S 3 1 (N) + 3S 2 (N)S 1 (N) 



+-S 3 (N) + (3"{N + 1) - A(3\N + 1)51 (iV) - 4/3 (JV + l)( 2 + Cs 



N3+8N2+UN+ 2 2 sm 



TV 4 - 2/V 3 + 5N 2 + 2N + 2 
2 ——_ . ___ — — -C2 



N(N + if (N + 2) 2 iy J (N — 1) N 2 (N + if (N + 2) 
7N 5 + 21N 4 + 13iV 3 + 21iV 2 + 187V + 16 

(N-1)N 2 (N + l) 2 (TV + 2) 2 ^ ^ 



5 



N 6 + 8N 5 + 23N 4 + 5AN 3 + 9AN 2 + 72N + 8 



-A 



N(N + if (N + 2) 3 
(N 2 - N - 4) 



Si (AO 



(N + 1) (N + 2) 



r/9'(iV+l) + 



P 3 (iV) 



(A 7 " — l)A r4 (A r + l) 4 (iV + 2) 4 



(36) 



PS, (2) 



(N) 



C F T, 



F 1 R 



— 4 



TV 4 + 27V 3 + 5iV 2 + AN + 4 
' (iV — l)iV 2 (iV + 1) 2 (A^ + 2) 
(N 2 + iV + 2) 2 



S 2 {N-l) 



c 2 + 



4 P 4 (iV) 



NS,(2) 

l gq,Q 



(N) 



C F T, 



F 1 R * 



(AT - 1)AT 2 (AT + 1) 2 (N + 2)" (JV - l)N 4 (N + 1) 4 (N + 2) 3 

/224 8 \ 40 8 

" (-27 + 3C2) S X {N - 1) + y S 2 (iV - 1) - -S 3 (iV - 1) 

2(3AT + 2)(AT- 1) 



•(37) 



<2 



3N(N + 1) 

(AT - 1)(219AT 5 + 428N 4 + 517N 3 + 512N 2 + 312N + 72) 

5AN 3 (N + l) 3 



(38) 



where 



P 2 (N) = 12A^ 8 + 54A^ 7 + 136A^ 6 + 218A^ 5 + 221A^ 4 + 110A^ 3 -3A^ 2 -24A^-4 (39) 
P 3 (N) = 2N 12 + 20N 11 + 8QN 10 + 192N 9 + 199N 8 - N 7 - 297N 6 -A95N 5 -51AN 4 

-488 N 3 - 416 N 2 - 176 N- 32 (40) 
P 4 (N) = N 10 + 8N 9 + 27 N 8 + 33N 7 - 71N 6 - 275 A^ 5 - 403 N 4 - 448 N 3 

-408N 2 - 208N - 48 . (41) 



Note that in the above expressions the equality 

Li 2 (x) 



M 



1 + x 



(N) - ( 2 (3(N) = (-1) 



N 



S_ 2A (N-l) + -( 3 



(42) 



can be applied. Therefore, the operator matrix elements depend on one non-trivial basic function 
[29,32] only and no sum with index { — 1} contributes. Concerning the non-trivial harmonic sums 
emerging the operator matrix elements are of the same complexity as the 2-loop anomalous 
dimensions. 

The longitudinal structure function Fl(x,Q 2 ) consists of the light and heavy flavor contri- 
butions 

F L (x,Q 2 ) = F^ ht (x,Q 2 )+F^(x,Q 2 ) 

= Cf s (x, a a , %r J <g> q N s(x, jj 2 ) + Cf \ x, a s , % J <g> q s (x, jj 2 ) 



g(x,fi 2 ) . 



(43) 



We choose Q 2 = fi 2 as uniform factorization scale. The Wilson coefficients in (43) are of the 
form 



x, a s ) + I x, a s , ^ 1 , i = S, NS, g . 



Q' 2 



(44) 



The heavy quark contributions are given by 



H L,g (%, a s, ^ 2 



rjPS / Q 
H L,q \ X , a s, ^ 



rrNS ( Q 2 
n L,g [ X , a s, m2 



-P^c { P ln(Q 

2 19 °L,q Ln 



771 



+ a 



11 99 



ln z 



5(2) 

(Q 2 



+n {2) + 77 (1) 
+a Qg + a Qg 



- p£ ] + m 



\my 2 

'L,q ^ ^ 19 



<3) 

2 / "i,g "T" °L,<? 



2-PS,(2) 3 
a s C L,q + °s 



_lp(0)p(0) ln 2(^_) 

+1 pPS,(i) i n f + ps>( 2 ) _ p(0) 
+ 2 <?<? m \ m 2j +a Qq a Qg^gq 



(1) ~PS,(3) 



-/3 0)QC «ln(^ j +c L> , 



~NS,(2) 



+ a: 



(45) 



(46) 



XC L,g + C L,q 



(47) 



The Wilson coefficients for heavy quark production consist of terms oc (m 2 / Q 2 ) k ,k > 0,/ceN and 
the logarithmic and constant contributions oc \n l (Q 2 /m 2 ), I > for on-shell massive quarks. The 
latter terms do not vanish in the limit m 2 — > and can be calculated solving the renormalization 
group equations for F® Q (x, Q 2 ). 



4 The Small- x Limit 



In the small x limit the heavy quark Wilson coefficient at 0(a s ), H^ g (x, Q 2 /m 2 ; fi 2 /m 2 ), vanishes 
with x since its leading pole is situated at N — — 1. One expects the following leading and next- 
to-leading small-x behaviour 



d {1) 

Ob 



oo 
k=2 



d 



(1) 



In 



* ~\x) {2) \n k - 2 (x) 
~r a k 



x 



+ 



X 



(48) 



In 0(a 2 ) the leading [24,33] small-x terms for ij 2 = Q 2 are 



45 



-S2QTR 



9 ' 



(49) 



with % — A,F for the gluonic and pure singlet contribution, respectively. As seen in Eqs. (8, 9, 
13, 14, 15), these terms stem from the small x behaviour of and Pfff, (21,22), and from 



c Lgi c Lg PS - The terms oc Tr for both contributions scale by the color factors, Ca, resp. Cp. Hl q 
contains an additional term at 0(a^), oc T| \n(Q 2 /m 2 



P^Lq' NS does not contain a term oc 1/x 



but is less singular for small values of x. The corresponding coefficients are obtained in finding 
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the contributions oc 1/(N — 1) in (45,46). Most of the respective functions were given above. 
We further note that c^ q (N = 1) = 2C F and c ( ^ q (N -> 1) oc -(32/9)C F T R N f /(N - 1), [9,27]. 
In 0(a 3 s ) the leading small a; contributions to H^'f^x, Q 2 /m 2 ; jj 2 /m 2 ) result from the small 



x terms in c Lq ^(x) only and are proportional to that of the light flavor contributions [13]. The 
remaining heavy flavor corrections are less singular. The leading terms cx \ jx are 

128 



d {2) 

a 2,A — 



-CaC{Tr 
6 

-32CaC f T r 
256 



34 / 



In 2 



Q 2 



10 



In 



28 
27 



27 
32 



C F T 2 R (2N F + 1) In 



m 2 



d (2) - 

a 2,F — 



+C F T 2 R {2N F 

-32C 2 F T R - 
F [3 

+32CaC f T r 
+C F T 2 R {2N F + 1) 



2756 65, 



64 

+ —CaC f T r 



56 



9 



- C2 - 4C 3 



+ C A T2(2iV F + 1) 




32 
~9~ 



101 



-8C2 



C2 - 4Cs 



Among the subleading terms not stemming from those oc T R scale by the color factor. 

H^ S q ^\x,Q 2 /m 2 ] fi 2 /m 2 ) is regular for N — 1. Similarly to the treatment in [34,35] one might 
consider its singular behaviour around N — 0, however, the small-x resummation for these terms 
were not yet derived. 



(50) 



(51) 



(52) 



5 Numerical Results 



In the following we will give some numerical illustrations of the effect of the heavy flavor con- 
tributions in the limit Q 2 ^> m 2 . These results are, unfortunately, of limited phenomenological 
use, since one expects, similar to the case of the NLO corrections [24], that these corrections 
become effective at large values of Q 2 ~ 1000 GeV 2 , where no data on F L are available at 
present. At lower values of Q 2 power corrections do still contribute. In the case of F 2 (x,Q 2 ), 
a sufficient description by the asymptotic expression could be obtained for scales Q 2 > 30 GeV 2 
already, cf. [24]. 

We illustrate the size of the contributions at NLO and NNLO for F®®{x,Q 2 ) choosing the 
parton distributions as follows at Ql = 30 GeV 2 



xq NS {x) 
xqps(x) 

xg(x) 



AWr aNS (l -x) 



Na as 



N g x a ° 



l-x) 

:i-x) 



+ c s x 



" 3 + CgX U 



(53) 
(54) 

(55) 



We apply this parameterization both for the NLO and NNLO effects for illustrative pur- 
poses. The parton distributions at higher scales are obtained by evolution. The values of 
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the parameters at Ql are listed in Table I. The strong coupling constant a s (Q 2 ) is cal- 
culated using the values of Aqcq NS,( ' 4 ' ) determined in [36]. One obtains: af LO (30 GeV 2 ) = 
0.1977, resp. 0.1708(Q 2 = 100 GeV 2 ), 0.1132(Q 2 = 10 4 GeV 2 ) and of NLO (30 GeV 2 ) = 
0.1928, 0.1673(Q 2 = 100 GeV 2 ), 0.1118(Q 2 = 10 4 GeV 2 ). The charm quark mass was chosen to 
be m c = 1.5 GeV. 



Parameter 


N t 


di 


h 


Ci 


di 


NS 


1.00000 


0.50000 


3.0000 






S 


0.60000 


-0.30000 


3.5000 


5.0000 


0.80000 


gluon 


0.11518 


-0.32230 


6.0445 


0.9618 


0.00422 



Table 1: The parameters of the quark non-singlet, singlet and gluon distribution at Q 2 = 30 GeV 2 . 

The asymptotic heavy flavor non-singlet contributions together with the light flavor terms 
are shown in Figure 1 and turn out to be small due to the shape of the input distribution and 
since they emerge only at 0(a 2 ). Their contribution shrinks with growing Q 2 . The asymptotic 
heavy flavor pure singlet part added to the light-flavor contributions of F L (x,Q 2 ) are depicted 
in Figure 2. Also here the the first contribution is obtained at 0(a 2 ), but the effect is much 
larger if compared to the non-singlet part due to the small-x behaviour of the corresponding 
distribution function oc x~ a3 . The contribution rises with Q 2 and amounts to ~ 0(1/5... 1/6) of 
the gluon contributions, shown in Figure 3. The gluon contribution is the largest and emerges 
already at LO. It grows towards the small x region. In all cases the NNLO result is larger than 
that at NLO. To obtain a complete picture for the heavy flavor contributions to Fl(x,Q 2 ) the 
corrections at scales of lower values of Q 2 have to be calculated. 



H 0.014 — 



O) 0.012 \r 
§ 0.01 \r 
kT" 0.008 F- 



Q 2 = 1 4 GeV 2 



Q 2 
Q 2 



1 2 GeV 
30 GeV 2 




X 



Figure 1: The light flavor and asymptotic heavy flavor non-singlet contributions due to charm to 
F L (x,Q 2 ) in NLO and NNLO. Upper lines: NNLO, lower lines NLO. 
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X 



Figure 2: The light flavor and asymptotic heavy flavor pure singlet contributions due to charm to 
F L (x,Q 2 ) in NLO and NNLO. Upper lines: NNLO, lower lines NLO. 




X 



Figure 3: The light flavor and asymptotic heavy flavor gluon contributions due to charm to F L (x, Q 2 ) 
in NLO and NNLO. Upper lines: NNLO, lower lines NLO. 
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6 Conclusions 



We have calculated the heavy flavor contributions to the structure function Fl(x,Q 2 ) in the 
region Q 2 ^> itlq at 0(a 3 s ). In this kinematic regime the respective terms are obtained as 
the logarithmic orders oc \n k (Q 2 /mq) and the constant term. Power corrections cannot be 
determined using the method of the present paper. In this approximation the heavy flavor 
Wilson coefficients are given by a convolution of the light-flavor Wilson coefficients and universal 
operator matrix elements, which contain the information on the heavy quarks. At NLO a 
numerical comparison of the complete calculation to the asymptotic case was possible and scales 
in the range Q 2 > 1000 GeV 2 were identified to apply the asymptotic relation for F®®{x,Q 2 ). 
This is likely to be the case at NNLO too. We presented numerical results for the asymptotic 
O(al) corrections added to the light-flavor contributions. The largest contribution is due to 
the gluonic term, followed by the pure singlet term, which is a factor ~ 5 smaller in the small 
x region. The flavor non-singlet contribution is very small. The leading small x terms of the 
heavy flavor Wilson coefficients H^ S q ^(x, Q 2 /m 2 ) were determined. In O(a^) the pure heavy 
flavor terms contribute to the next-to-leading small x terms. In part of the terms scaling by 
Cf{Ca) is observed comparing the respective gluonic and quarkonic contributions. 
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7 Appendix 



The Mellin transforms used in the present calculation may be found in [27]. Some expressions 
can be written in the more compact form given below. 



M[ln(l + z)](A0 = l{ln(2)-(-l) JV [5_ 1 (iV) + ln(2)]} 



M[hi(z) ln(l + z)}(N) 
M[ln 2 (^) ln(l + z)](N) 



M[U 2 (-z)](N) 

M[\n(z)U 2 (-z)](N) 
M[U 2 (-z) + \n(z) ln(l + z)](N) 
M[U 3 (-z)](N) 
M[$i(z)](N) 



N 

l[ln(2)-/3(7V + l)] 
-^N2)-/3(iV + l)]-i/3'(7V + l) 



(56) 
(57) 



AT3 



{ln(2)-(-l) JV [5_ 1 (iV) + ln(2)]} 



■1) 



A? 



N 2 



S- 2 (N) + 



<2 



"I) 



A? 



N 



A [i n(2 ) - (3{N + 1)] + JLp(N + 1) - ±(3"(N + 1) 



J + -1 [ln(2) - 0(JV + 1)] 



2AT A^ 2 

A-lpn(2)-^ + l)]-i^ + l) 
[Ca + 2F(N + 1)] 



2AT 

^C3 + ^C 2 -^[ln(2)-/3(7V + l)] 
-1) 



(58) 

(59) 
(60) 
(61) 
(62) 



Af+l 



iV 



{25' 1) _ 2 (iV)+C2[5 , i(JV)-S'_ 1 (JV)]} 



■1) 



N+l 



N 



C 2 ln(2) 



(63) 



1 

iV 



|2M 


Li 2 (x) 




_ 1 + x _ 



W-^C 2 + ^iW+3C 2 /3(iV + l) 



Cs 



+2S 1 (iV)/3 , (iV + 1) - /3"(jV + 1) +^ - C 2 ln(2) 



where 



$i(z) = 2Li 2 (-z) ln(l + z) + ln 2 (l + s) ln(z) + 2S 1)2 (-z) . 
M [Li 2 (x)/(1 + x)] (N) is a basic function, cf. [28,30]. 



(64) 



(65) 
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